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We show that the entanglement spectrum associated with a certain class of strongly correlated 
many-body states — the wave functions proposed by Laughlin and Jain to describe the fractional 
quantum Hall effect — can be very well described in terms of a simple model of non-interacting (or 
weakly interacting) composite fermions. 


Understanding of strongly correlated many-body 
quantum states has been significantly enhanced in recent 
years by the introduction of the entanglement spectrum 
(ES), i.e. the spectrum of eigenvalues of the reduced 
density matrix of a subsystem. 1 The scrutiny of the ES 
has been particularly fruitful in the study of fractional 
quantum Hall (FQH) wave functions — the archetype for 
topological phases of matter resulting from strong inter¬ 
action effects — where this analysis often reveals a special 
underlying entanglement structure useful for identifying 
and classifying different topological phases. 1-6 

For FQH states in particular it has been established 
that there is a deep connection between the entanglement 
structure of bulk wave functions (as seen in their ES 7-9 ) 
and the structure of edge excitations. 10-13 Central to this 
connection so far have been arguments based on the com¬ 
mon conformal held theory description of the bulk and 
edge physics, which is known to apply for the class of 
wave functions that can be constructed in terms of ex¬ 
pansion functions known as “conformal blocks”. The ES 
for this class of what we shall refer to as “simple” FQH 
states has been extensively studied using the machinery 
of conformal held theory. 4,5,9 

Largely separate from these developments, it has been 
well documented that the composite fermion model of 
the FQH effect can successfully account for many exper¬ 
imentally observed features of the effect, particularly in 
the lowest Landau level. 14 It does so by positing that 
they are due to the integer quantum Hall (IQH) effect 
of non-interacting (or weakly interacting) quasiparticles 
known as composite fermions. While composite fermion 
theory presents an appealingly simple physical picture to 
explain the FQH effect, it has proved to be very chal¬ 
lenging to relate the many-body composite fermion wave 
functions to constructs in conformal held theory. While 
such constructs do exist, 15-18 they are somewhat com¬ 
plicated, and this has so far limited the effectiveness of 
conformal held theory machinery in describing the ES 
of the Jain states. In this sense the Jain states are not 
“simple” FQH states. 

In this work we present a very different approach to 
constructing entanglement spectra that bypasses the po¬ 


tential difficulty associated with explicitly writing wave 
functions in terms of conformal blocks. We show that 
the low-lying (highest weight) part of the ES of the Jain 
FQH states can be accurately described by a modified 
ES of non-interacting (or weakly interacting) composite 
fermions in filled Landau levels (in other words a mod¬ 
ified ES of the IQH states 19-21 ). To demonstrate the 
effectiveness of our method, we shall focus on presenting 
results for two fundamental examples, namely the Laugh- 
lin state at filling factor v = 1/2 for bosons (a simple 
FQH state that was previously studied using conformal 
field theory methods 9 ) and the Jain state at v = 2/3 for 
bosons (a non-simple FQH state). (The method applies 
equally well to fermionic FQH states.) Our results for the 
Laughlin case are in good agreement with the approach 
based on conformal field theory, but the real advantage 
of our method becomes evident for the Jain states, where 
the conformal Held theory approach appears much more 
complicated, and has therefore not been worked out. 

In the spirit of the original composite fermion model, 14 
we start with the ES for the IQH states, 20,21 treating it 
in a general framework for determining the ES of non¬ 
interacting systems. 19 We shall briefly review its deriva¬ 
tion here. For simplicity we shall consider spinless par¬ 
ticles. Also, to remove any additional complications due 
to edge physics we consider a system without physical 
edges. A standard technique to achieve this is to solve 
the problem on the surface of sphere. 22 Here we have an 
integer n Landau level problem (where the filling fac¬ 
tor is v = n) for N particles on a sphere of radius 
VQ (in units of magnetic length) that encloses a ficti¬ 
tious magnetic monopole of strength 2 Q = (N — n 2 )/n. 
There is rotational symmetry about the z-axis, which 
leads to the single-particle orbitals of the problem be¬ 
ing labelled by the z-component of angular momentum 
in = — (2Q + n)/2, — (2 Q + fj)/2 + 1,..., (2 Q + cr)/2 in 
addition to a Landau level index a = 0,1,..., n — 1. The 
single-particle orbitals are 4>m,a{ r), where r lies on the 
surface of the sphere (technically, (f> m ,a( r) are monopole 
harmonics. 23 The monopole harmonics also depend ex¬ 
plicitly on Q , but for simplicity we suppress this depen¬ 
dence in our notation). 
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In order to connect to existing calculations we shall 
focus on the entanglement spectrum for wave functions in 
real space (as opposed to e.g. momentum space), which 
leads to the so-called real-space entanglement spectrum 
(RSES). 4 - 6 

We shall consider cuts along lines of latitude, so that 
that rotational invariance about the z-axis is preserved 
and the ^-component of angular momentum L z remains 
a good quantum number. Once the system is cut, L z is 
bipartitioned as L z = Lf + Lf. Additionally, the total 
particle number is bipartitioned as N = IV )4 + Nb. 

For the full fermionic iV-particle system, we have a 
basis for the Hilbert space that consists of Slater deter¬ 
minants built from the orbitals 4>m,a- For the A and B 
subsystems, the restrictions of these orbitals to the re¬ 
gions A and B still form a complete single-particle basis, 
although the restricted orbitals must be renormalized to 
account for the fact that only part of the single-particle 
weight is in each subsystem. Once restricted to a sub¬ 
system, the orbitals ^ mij(Tl and 4>m 2 ,a- 2 are orthogonal 
whenever mi ^ m 2 (because the cut is made to preserve 
L z ), but are generally not orthogonal whenever cy 02 
(though they remain linearly independent in this case). 
The Fock space of the A and B subspaces is spanned by 
Slater determinants in terms of these restricted orbitals. 

We can now write the Schmidt decomposition of the 
state IV’) as 


IV’) = J2 e 5,72 K A > 0 K B ) C 1 ) 


where \ipf) and |V^) belong to the Fock spaces for the A 
and B subsystems that we have just described. The non¬ 
zero Schmidt coefficients, e - ^/ 2 , are conveniently writ¬ 
ten in terms of the entanglement energies, Vi- The index 
i labels the vectors in the Schmidt basis for the A sys¬ 
tem (at least those with non-zero Schmidt coefficients). 
Since Na and Lf are conserved by the cut, we can choose 
the Schmidt basis to consist of eigenstates of Na and 
Lf. We can then write Va r A ,L A ,i f° r the entanglement 
energy of the i th Schmidt state with given Na and Lf. 
In fact, with any Schmidt state with Na particles and 
given Lf, we can associate an AA-tuple of single par¬ 
ticle momenta m = [m\,... ,m^ A ) with Lf = fZp m p 
and m± < m 2 < ... < m^ A . If there is only a sin¬ 
gle Landau level involved, the Schmidt state is just the 
unique Slater determinant in the Fock space of system A 
labelled by m and we can replace the label i above by 
m. With multiple Landau levels, there will be a number 
of Schmidt states associated with each m. These states 
are superpositions of the Slater determinants with this m 
and different choices of the numbers of particles in each 
Landau level. 

One can find exact expressions for the entanglement 
energies and Schmidt states in the case of a non¬ 
interacting Landau level problem, 20 , using an argument 
proposed by Peschel. 19 If one constructs the correlation 


matrix, 


Cm 1 ,/7 1 ,m 2 ,o-2 / V > mi,(Ti ( r )^ , m 2 ,o- 2 ( r ) dr 

J A 

= dm 1 ,m,2 / 4 > m 1 ,a 1 (. r )4 > m 1 ,a 2 { r ) dv (2) 

J A 

then the eigenvalues, A mCT , of C are related to the entan¬ 
glement energies by 


& = Z)' 


'm,cr,i c m,cr 


constant, 


( 3 ) 


where the single-particle entanglement energy function , 
Cm.a, is defined by 


Cm, a — ffig 


1 Am, 

Am,cr 


( 4 ) 


Note that now a labels the eigenstates of the correlation 
matrix at given m, rather than a Landau level. The full 
Schmidt states are Slater determinants built from these 
eigenstates and o mjOVi denotes the occupation number of 
the correlation matrix eigenstates labelled by (m, u) in 
the Schmidt state with label i. Examples of the RSES 
for the v — 1 and v = 2 IQH states calculated using 
this method are included in Fig. la and Fig. 2a for later 
comparison with our results for the FQH case. 

We now come to the statement of our main result: 
For the class of FQH states proposed by Laughlin and 
Jain to describe the FQH effect in the lowest Landau 
level at filling factors v = n'/ (2 n' + 1) [for fermions] and 
v = n!/ (n r + 1) [for bosons] with n' integer, the low-lying 
V in the associated RSES can be accurately described by 
the same model as the IQH states (Eq. 3), but a differ¬ 
ent single-particle entanglement energy function (Eq. 4). 
In addition, when multiple effective Landau levels are 
present in the Jain wave functions, one also needs to 
take into account a simple exchange-entanglement-energy 
term that we shall shortly describe with an example. 

In this connection the Landau level index a for the IQH 
wave functions becomes the effective Landau level index 
of the Jain states, a' = 0,1,..., nl — 1, the magnetic 
field described by the monopole strength Q is replaced 
with an effective magnetic field and an effective monopole 
strength Q' = (N — n' 2 )/n' and the occupation basis of 
single-particle Landau level orbitals now becomes an oc¬ 
cupation basis of so-called composite fermion orbitals in 
the effective Landau levels, whose single-particle eigen¬ 
states are labelled by ml = —{2Q'-\-a')l2,—(2Q'+(j')/2-\- 
l,...,(2Q' + a')/2. 

It is useful to expand e TOjCr from Eq. 4 as a power series 
in m about the mid-point m = 0, as it turns out that 
we need to keep only very few terms in this expansion 
for a very accurate description of e mjCT . For example, in 
the case of a single Landau level (u = n = 1 for spinless 
particles) or a single effective Landau level {n' = 1 and 
replacing m with ml in the expansion) 

Cm,a ~ ao + aim + a 2 m 2 + a 3 m 3 + ..., (5) 
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FIG. 1. (color online), (a). Exact RSES for the v = 1 
state for N = 50 particles, with equal-size A and B regions 
and Na = 25. Embedded number labels indicate the degen¬ 
eracy of £ (up to ALf = 8). 24 (b). Fit of the single-particle 
entanglement energy model (Eq. 5 in terms of the composite 
fermion angular momentum m' and truncated at order m' 3 ) 
to the ES of the v = 1/2 Laughlin state calculated numeri¬ 
cally for the same cut and system size using the method in 
Ref. 25. Fitting parameter data is included in Appendix B. 
Entanglement energies £ and angular momentum labels Lf 
are relative to the lowest lying state £o in the lowest angular 
momentum sector Lq i.e. A£ = £ — £o and A Lf = Lq — Lf. 


where, for the IQH case, do, ai etc. are functions that 
can be expressed in terms of the eigenvalues of the matrix 
C and its derivatives (note that they depend on the sys¬ 
tem size, via their dependence on the monopole strength 
Q). For 2 filled Landau levels (v = n = 2 for spinless 
particles) or two effective Landau levels ( n' = 2) the ex¬ 
pansion becomes instead 

e m ,cr ~ a 0 ,o + a 0 ,im + a 0 ,2 m 2 + a 0y3 m 3 + ... (6) 

+ [<j — 1/2] {ai j0 + aipm + ai |2 rn 2 + ...}. 

In our approach for the FQH case the coefficients (ij t k 
now play the role of fitting parameters in a model given 
by truncating the single particle energy function at a low 
degree in rn!. Their values are determined by fitting this 
model to the RSES generated numerically from a micro¬ 


scopic construction of the FQH state using other meth¬ 
ods. Further details of the fitting procedure are given in 
Appendix A. To obtain an appropriate numerical spec¬ 
trum to fit to for large system sizes we use the method 
proposed in Ref. 25 (a method that works more efficiently 
for bosonic rather than fermionic states). 

To illustrate how our approach works we present two 
examples of applications to bosonic FQH states, where 
we choose symmetrical, equal-sized regions A and B (i.e. 
a cut along the equator of the sphere) and Na = Nb = 
N/ 2 . 24 The method also applies for more general cuts. 

As our first example, we consider at the RSES of the 
Laughlin state of bosons at filling factor v = 1/2 (this 
state also occurs in the Jain series for n' = 1). The 
RSES for v = 1/2 is accurately described by the single¬ 
particle energy function in Eq. 5, written in terms of 
the composite fermion angular momentum labels m! and 
truncating the expansion at order ?n ,3 . In this truncated 
expansion, the coefficients a±, a 2 and <23 are free fitting 
parameters. It is not necessary to include a constant 
term (ao) in this fit since its value can be fixed by nor¬ 
malizing the spectrum (see Appendix A). An example fit 
to the numerically-calculated RSES of the v = 1/2 state 
is shown in Fig. lb. Fitting parameter data is included 
in Appendix B. We find that by far the most dominant 
contribution comes from the aq term, indicating that the 
spectrum is almost linear. In Fig. la we also show the 
RSES of the IQH state for v = 1 for comparison. 

The Laughlin state falls into the class of simple FQH 
states and its RSES has been studied previously using 
conformal field theory techniques. In particular Ref. 9 
describes in detail a method to model the Laughlin RSES 
in terms of the energy levels of an entanglement Hamil¬ 
tonian that is given by writing down every allowable con¬ 
formal field theory operator (in this case chiral boson op¬ 
erators, their derivatives and powers) order-by-order in 
their scaling dimension (with terms at lower scaling di¬ 
mension providing more relevant contributions). These 
operators come with unknown coefficients that can be 
fitted to the numerically-calculated RSES using a similar 
procedure to that described in Appendix A. In Ref. 9 it 
was shown that the RSES of the v = 1/2 state can be very 
accurately described by truncating at scaling dimension 
2, leading to an entanglement Hamiltonian containing 3 
terms (and therefore 3 unknown fitting parameters). We 
have checked that for the same data set as used here 
(Fig. lb) the fit of this truncated entanglement Hamil¬ 
tonian to the RSES is in excellent agreement with our 
non-interacting model (in fact, the quality of the fit as 
defined in Appendix A is typically improved by an or¬ 
der of magnitude compared to our approach). One could 
in principle reproduce our model directly (with the in¬ 
clusion of additional interaction terms) by fermionizing 
Dubail-Read-Rezayi’s entanglement Hamiltonian. The 
fact that our result agrees closely with the numerically- 
calculated RSES implies that any additional interaction 
terms must only be small corrections, which justifies the 
assumptions made by the composite fermion approach. 
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FIG. 2. (color online), (a). Exact RSES for the v — 2 state 
for IV = 48 particles, with equal-size A and B regions and 
Na = 24. (b). Fit of the single-particle entanglement energy 
model (Eq. 8 ) to the v = 2/3 Jain state calculated numerically 
for the same cut and system size using the method in Ref. 25. 
Fitting parameter data is included in Appendix B. <j' tot labels 
the number of particles in the upper effective Landau level 
for each /. Entanglement energies / and angular momentum 
labels L^ are relative to the lowest lying state £o in the lowest 
angular momentum sector Lq i.e. A£ = / — /o and A Lf = 
Lq — Lf. Number labels embedded in the plots indicate the 
degeneracy of £ where ambiguous . 24 


Our second example is the RSES of the Jain state of 
bosons at filling factor v = 2/3 (this state occurs in the 
Jain series for n! = 2 for bosons). Our technique was 
not able to improve upon conformal field theory in treat¬ 
ing simple FQH wave functions, but where it really ex¬ 
cels is in its treatment of the Jain states. In this case 
we now use the single-particle energy from the v = 2 
IQH state (Eq. 6 ), but we augment it to take into ac¬ 
count an exchange-likc interaction (treated at the level 
of “mean-field theory”) for fermions in different Landau 
levels. This term can alternatively be thought of as a 
“charging energy”. Our ansatz to describe the RSES of 


the v — 2/3 Jain state is 

^NA,Lf,i = ^ T C (ATV) , (7) 

m,<j 

with the single-particle energy function given by 

e m ', a' ~ a 0 ^m' + [o' - 1/2] {aqo + aiqffi'} , (8) 

and the exchange-interaction, or charging energy term, 
given by 


AN= J2 n m ,^[o'-1/2}. (9) 

m' ,a' 

This ansatz requires 4 fitting parameters ao,i> ai, 0 ) 
and c (once again, a constant term a 0,0 is not fitted be¬ 
cause it can be fixed by normalization). 

In Fig. 2b we show an example fit to the RSES of the 
2/3 state for a large system of N = 48 particles. Fitting 
parameter data is included in Appendix B. We observe 
that the most dominant contribution arises from the coef¬ 
ficient ai, 0 j which can be thought of as an entanglement 
“cyclotron energy” term. We also find that the coeffi¬ 
cient of the exchange interaction, c, is positive (so due 
to this additional term, the branches in the RSES for 
1 / = 2/3 are further apart than they would have been 
otherwise). Note that there are degeneracies in the fit¬ 
ted entanglement energy eigenvalues, but these can be 
lifted by truncating the entanglement energy expansion 
at higher order in m !. For comparison, in Fig. 2a we plot 
the RSES of the v = 2 IQH state. 

To summarize, we have described how a RSES of non¬ 
interacting composite fermions can be constructed to ac¬ 
curately approximate the RSES for certain strongly cor¬ 
related FQH states. Key to this construction is the ob¬ 
servation that the single-particle entanglement energy 
function underlying the description of the RSES for the 
IQH states can be simply modified in order to describe 
the RSES of FQH states, treating the many-body RSES 
within a non-interacting approximation and allowing for 
very basic mean-field theory like exchange energy correc¬ 
tions in the multi-Landau level case. We note that this 
description of the RSES closely parallels the description 
of the real energy spectrum of the edge of a Hall droplet 
in terms of composite fermions (see e.g. Ref. 26). 

The quality of this approximation could be improved 
by allowing for higher order corrections to the underly¬ 
ing entanglement energy function (at the cost of need¬ 
ing more fitting parameters). These corrections can 
be thought of as allowing for additional inter-composite 
fermion interactions at the level of “mean field theory”. 
Alternatively, weak interaction corrections can be added 
to the model perturbatively, working in the framework of 
degenerate perturbation theory as applied to the entan¬ 
glement Hamiltonian. However, we find that such cor¬ 
rections only provide a marginal improvement to the fit. 
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Appendix A: Fitting algorithm 

In this appendix we shall briefly describe our procedure 
for evaluating the goodness of fit of the RSES to the 
single-particle models proposed in this work. 

The problem is a many-parameter optimization where 
we aim to minimize a weighted sum of squared differ¬ 
ences between corresponding £ in the lowest lying part of 
the RSES and single-particle energy spectrum for a given 
set of fitting parameters {a,j,k}- By corresponding £ we 
mean that, for comparison, we normalize both the model 
spectrum and the numerically calculated spectrum such 
that the lowest lying £ in the lowest A Lf sector is set to 
zero. Then we order the set of £ in the model spectrum 
(call them £ model ) and the lowest lying part of the RSES 
(call them £ RSES ) by their A Lf values, then for each sec¬ 
tor we sort in order of increasing £ RSES or £ model value and 
finally we take the sum of squared difference between the 
lowest £ rses and lowest £ modt>1 , the difference between the 
next lowest £ RSES and next lowest £ mode i and so on. This 
sum can also be weighted in various ways, for instance 
if we want to give increased importance to matching up 
states with the lowest values of £. In general, therefore, 
we aim to minimize a fitting function of the form 

R( K/J) = E [&,.*> - 2 w it (Ai) 

i 

where {a,j,k} denotes the set of free parameters in the 
model, and the index i appearing in the sum denotes the 
ith value of £ model , £ RSES within the ordered set (ordered 
in the sense described above). The factor W, assigns an 
optional weight to each term in the sum. In addition, we 
only include entanglement spectrum eigenvalues below a 
certain specified cut-off sector AL^ t _ off , i.e. Wi = 0 if it 
refers a state with A> A L^ ut off . 

In our application we always include a factor in Wi that 
divides out the number of states in each sector (i.e. for 
every £* H labelled by the same A Lf we divide by a factor 
of N al a that counts the total number of £ model with that 
A Lf). The reason for doing this is because otherwise 
the fitting function would assign overwhelmingly more 
weight to fitting the higher sectors (because the counting 
of states in the RSES grows superpolynomially). Once 
the N al a factor is divided out, each sector counts for the 
same total weight in Eq. Al, i.e. Y^i„ahAL A W* = 1 for 
all A Lf. Another aspect to consider is that, in the RSES, 
linearly smaller £ correspond to exponentially greater co¬ 
efficients in the Schmidt decomposition. In this sense, it 
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should be more important physically to match up the 
lower lying £. For this reason, we might also consider 
including a factor exp(—£^ SES ) in our weight function. 
This addition works well for the Laughlin case, however 
the Jain case is more complicated due to the presence 
of the branch structure in the RSES. Consequently, for 
studying the Laughlin state we use the fitting function 


E 

i\ AL?<AL a . « 

1 z — cut-off 


(q SE s-Co d .i y 

Na if 


whereas for fitting the Jain state we find that for fitting 
all branches with equal weight it is better to use 


R,„r, = 


E 

i\AL?<AL? ut _ oit 


(&SES ~ 

Nal* 


The quality of a given fitting model (e.g. a given set of 
fitting parameters {ay*,}) can be assessed by the minimal 
value of the fitting function obtained for that model. The 
lower the minimal value of the fitting function, the better 
the quality of the fit. 

In order to solve the minimization problem we use the 
Powell method provided by Scientific Python (SciPy) ver¬ 
sion 0.11 and above, which involves a sequential ID mini¬ 
mization of each fitting parameter. The Powell method is 
found to be numerically stable for this problem. One also 
has to take steps to avoid finding the local rather than 


the global minimum, which we achieve by running the 
procedure a large number of times with different random 
starting parameters. 

Appendix B: Fitting parameter data 

In this appendix we tabulate the values of the fitting 
parameters cij^ obtained in the single-particle energy fits 
shown in Fig. 1 and Fig. 2. We also include the corre¬ 
sponding values of the fitting function R,{{a : j^}) defined 
in Appendix A to assess the relative quality of those fits. 

Example fitting parameter values and i? La u g hii„ 
(Eq. A2) values in the single-particle model proposed in 
this work for the Laughlin v = 1/2 state. Correspond¬ 
ing results plotted in Fig. 1. Smaller values of i? L aug hii„ 
indicate a better quality of fit: 

O-l _02_ 03 R L au g hlin 

0.243 0.00130 0.00214 8 8.12 x 

10" 5 

Example fitting parameter values and i? Tain (Eq. A3) 
values in the single-particle model proposed in this work 
for the Jain v = 2/3 state. Corresponding results plotted 
in Fig. 2. Smaller values of i? Jain indicate a better quality 
of fit: 

Qq,i 01,0 01,1 _ c AL^ t _ off i? Jain 

0.9279 3.371 0.6429 0.1557 8 0.945 



